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摘 要 : 基于 章 古 台地 区 22 块 樟 子 松 ( Pinus sylvestris var. mongolica) 人 工 纯 林 标 准 地 的 702 棵 樟 子 松 立 木 数据 , 构 


建 了 樟 子 松 固沙 林 冠 幅 一 胸径 关系 的 基础 模型 .广义 模型 及 基于 混合 效应 的 基础 模型 和 广义 模型 ;比较 了 随机 选 
择 样 本 木 .选择 平均 胸径 树 `. 选 胸径 较 小 树 和 选 胸径 较 大 树 4 种 方案 ,计算 混合 模型 随机 参数 时 的 混合 模型 预测 精 
度 ;最 后 分 析 了 不 同 林木 因子 和 林 分 变量 对 冠 幅 一 胸径 关系 的 影响 。 模 型 评价 指标 包括 决定 系数 ( 尼 ) .平均 绝对 
误差 ( MAE) 以 及 均 方 根 误 益 (RMSE ) 。 结 果 表 明 : 枝 下 高 (HCB) 、 相 对 植 距 ( RS) 和 林 龄 (4) 对 冠 幅 一 胸径 关系 影 
响 最 为 显著 ;混合 模型 拟 合 精度 (基础 混合 模型 R MAE 和 RMSE 分 别 是 :0.703 0.0.386 6 和 0.515 4; 广 义 混 合 模 


#0 R? MAE #1] RMSE 为 :0.705 1.0.382 2 和 0.513 6) 高 于 最 小 二 乘法 回归 (OLS) 模 型 (基础 模型 RR、MAE 和 RMSE 


分 别 为 :0.587 5 0. 469 6 0. 607 5; XA! R? MAE 和 RMSE 分 别 为 :0.661 8 0.415 5 和 0.550 0) 。 基 础 混合 模 
型 和 广义 混合 模型 差异 较 小 (2 模型 RR MAE 和 RMSE 均 相 差 1% 左右 ) 。 冠 幅 随 HCB 和 4 的 增 大 而 减 小 , 随 RS 
的 增 大 而 增 大 。 进 行 冠 幅 预测 时 ,推荐 使 用 基础 混合 模型 并 从 每 块 标准 地 选择 2 棵 平均 木 冠 幅 计 算 其 随机 参数 ， 


或 使 用 方法 较为 简单 的 OLS 广义 模型 预测 单 木 冠 幅 大 小 。 


关键 词 : 榜 子 松 ; 人 工 固 沙 林 ; 冠 幅 一 胸径 ; 混合 效应 ; 基础 模型 ; 广义 模型 


冠 幅 是 预测 树冠 表面 积 .树木 生物 量 、 林 分 郁 闭 
及 树冠 形状 的 重要 变量 " ,是 林 分 模型 的 重要 参 
数 ” ,可 用 于 计算 某 些 植被 的 防风 固沙 效应 指标 
(ASE TLR) ,其 结构 是 种 群 地 上 部 分 生命 活动 最 


差异 和 林 分 效应 的 广义 模型 进行 预测 。 
近年 来 ,混合 效应 模型 在 处 理 重复 测量 数据 和 

纵向 数据 方面 的 优势 被 不 断 发 所 ,该 类 模型 即 可 

描述 冠 幅 和 胸径 关系 在 不 同 林 分 间 的 差异 ,还 具有 


活跃 .影响 环境 最 激烈 的 部 分 。 定 量 认 知 树木 冠 
幅 对 林 分 的 密度 控制 .生态 效益 和 生产 能 力 的 最 优 
化 管理 具有 重要 意义 。 但 是 , 冠 幅 的 调查 成 本 较 
高 且 费 时 费力 ,于 是 建立 能 够 准确 预测 冠 幅 生长 变 
化 的 模型 便 十 分 必要 。 

目前 常用 的 冠 幅 预测 模型 中 ,一 种 是 只 含有 胸 
径 1 个 自 变量 的 简单 模型 , 男 一 种 是 由 胸径 和 其 他 
能 够 反映 林木 和 林 分 间 差 异 的 林木 因子 和 林 分 变量 
构成 的 广义 模型 。 有 学 者 认为 , 仅 含 胸径 1 个 预测 
变量 的 模型 既 可 以 满足 需要 ,还 具有 预测 变量 少 的 
应 用 优势 ““” 。 也 有 研究 表明 : 冠 幅 大 小 受 林 分 密 
度 、 立 地 条 件 、 和 森林 类 型 和 管理 措施 等 因素 的 影 
响 " ,简单 的 冠 幅 一 胸径 回归 模型 在 用 于 总 体 平 
均 水 平 预测 时 可 能 效果 较 好 ,但 应 用 于 个 体 预 测 时 ， 
在 大 密度 林 分 中 会 出 现 预 测 值 偏 大 ,在 稀 琉 林 分 中 
预测 值 偏 小 的 问题 '” ,因而 , 需 使 用 能 够 表示 林木 
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变量 选择 自由 以 及 拟 合 与 预测 精度 高 等 优点 ,在 林 
业 建 模 领 域 应 用 广泛 "目前 国内 外 已 有 一 些 
关于 混合 效应 冠 幅 预测 模型 的 报道 ,但 是 各 树种 间 
存在 差异 。Sharma “°°” 以 欧洲 云 杉 (Picea 
abies) 欧洲 山 毛 样 (Fagus sylvatica ) 和 欧洲 赤松 
(Pinus sylvestris ) 为 研究 对 象 ,建立 了 样 地 水 平 混合 
模型 ,确定 了 3 个 树种 冠 幅 模型 中 的 树木 因子 和 林 
分 变量 ,并 分 析 了 所 加 变量 对 冠 幅 预测 的 影响 ;也 有 
学 者 通过 使 用 2 水 平 (立地 与 样 地 ) 非 线性 混合 模 
型 ,建立 了 杉木 ( Cunninghamia lanceolata ) 单 木 冠 幅 
预测 模型 2 ; Fu EO 以 山西 地 区 的 华北 落叶 松 
(Larix principis-rupprechtii ) AU RMA, RIGA WE 
变量 的 2 水平 (地 域 与 样 地 ) 非 线性 混合 模型 可 以 
对 不 同 地 域 和 样 地 中 的 华北 落叶 松 冠 幅 进行 预测 。 
总 体 看 来 ,国内 利用 混合 模型 预测 冠 幅 的 报道 还 比 
较 少 , 且 多 为 混合 模型 与 传统 模型 预测 精度 的 比较 ， 
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关于 此 类 混合 模型 在 不 同 抽样 方案 条 件 下 预测 精度 
比较 的 报道 还 较 少 。 

樟 子 松 ( Pinus sylvestris var. mongolica ) 具有 而 
早 耐寒 、 耐 贫 盗 、 抠 干 通 直 、 生 长 快 . 防 风 固 沙 能 
强 等 优点 。 自 20 世纪 50 年代 被 成 功 引种 到 辽宁 省 
章 古 台 以 来 , 樟 子 松 在 中 国 “ 三 北 " 地 区 得 到 了 大 面 
积 的 推广 , 现 已 成 为 中 国 北方 干旱 半 干 旱地 区 造林 
的 主要 树种 "7" 。 据 研究 ,由 于 水 热 条 件 的 差异 , 辽 
宁 省 章 古 台地 区 的 沙 地 樟 子 松 与 天 然 分 布 区 ( 呼 伦 
贝尔 沙 地 ) 的 沙 地 樟 子 松 在 林 分 结构 和 生长 特性 方 
面 存在 较 大 差异 "" ,因此 ,在 经 营 管理 方面 也 不 同 
于 天 然 分 布 区 。 然 而 , 沙 地 樟 子 松 经 营 管理 方面 的 
文献 非常 萎 乏 , 据 笔 者 了 解 ,目前 榜 子 松 人 工 林 冠 
幅 一 胸径 模型 以 及 基于 混合 效应 模型 两 者 关系 的 报 
道 还 寥寥 无 几 。 鉴 于 此 ,本 文 以 章 古 台 沙 地 榜 子 松 
人 工 林 为 研究 对 象 , 以 期 建立 适合 该 地 区 樟 子 松 的 
冠 幅 预测 模型 ,确定 使 用 混合 模型 预测 冠 幅 时 最 优 
抽样 方案 ,为 沙 地 樟 子 松 人 工 林 的 合理 经 营 提供 理 
论 基 础 。 


1 材料 与 方法 


1.1 研究 区 概况 

研究 区 位 于 科尔沁 沙 地 东南 缘 的 辽宁 省 彰 武 县 
章 古 台地 区 (42° 39’ ~ 42。43' N, 122° 23' ~ 
122°33'E) ,海拔 高 度 为 (225 +5) m。 该 地 区 地 处 
中 温带 ,大 陆 性 季风 气候 ,年 均 降 水 量 500 mm 左 
右 , 年 蒸发 量 平均 在 1 600 mm 左右 。 土 壤 主要 类 型 
为 生 草 风沙 土 和 流动 风沙 土 。 具 有 流动 沙丘 、 半 固 
定 沙丘 PARED Er AF Be HS He SS, FR HE 
HEAL PAB Ze BK (Potentilla anserina ) P E fà F 
草 ( Cleistogenes chinensis ) , /|\ H $ ( Artemisia frigi- 
da) $ ty Bt (Arundinella anomala ) il] Ù ( Armeniaca 
sibirica ) 兴安 胡 枝 子 (Zespedeza daurica ) hE ( Arte- 
misia halodendron ) 、 大 果 榆 ( Ulmus macrocarpa ) 、 黄 
fil (Salix gordejevii ) PIPE ZE RE ( Agriophyllum squarro- 
sum), 研究 区 营建 模子 松 人 工 林 前 为 固定 
半 固 定 沙 地 。 
1.2 数据 来 源 

在 充分 踏 查 基 础 上 选择 合适 的 标准 地 ,标准 地 
的 选择 要 满足 3 个 原则 :首先 ,要 尽 可 能 覆盖 榜 子 松 
固沙 林 不 同 的 林 分 年 龄 (由 于 幼林 林木 冠 幅 较 小 ， 
容易 测量 ,因此 本 研究 在 建立 模型 时 未 包含 幼 龄 


林 ) ;其 次 ,要 覆盖 该 地 区 不 同 的 立地 条 件 ;第 三 , 标 
准 地 间 的 距离 不 少 于 50 m ,避免 紧 靠 道 边 .农田 , 排 
除 混交 林 。22 块 樟 子 松 人 工 纯 林 临时 标准 地 满足 
以 上 条 件 ,每 块 标准 地 的 规格 为 20 m x20 m。 标 准 
地 的 林 分 年 龄 涵盖 20 ~ 60 a 年 龄 段 ,包括 了 樟 子 松 
中 龄 林 及 以 后 的 各 个 阶段 。 立 地 覆盖 了 平坦 沙 
地 及 各 个 坡 向 的 坡 上 . 坡 中 和 坡 下 。 于 2016 年 4 月 
对 所 选 标准 地 内 的 立木 进行 了 每 木 检 尺 ,测量 过 程 
中 ,使 用 精度 为 0. 1 m 的 红外 线 测 树 仪 ( 快 特 能 
RD1000 ) W E (H) 和 第 一 活 枝 下 高 (HCB) ;使 
用 胸径 尺 测量 树木 的 带 皮 胸 高 围 长 ,用 以 计算 胸径 
( DBH) ;使 用 卷 尺 测量 树木 的 东西 冠 幅 (Cav ) 和 
南北 冠 幅 (CW\s ) , 取 两 者 的 平均 值 作 为 树木 冠 幅 
(CW) 值 ;通过 查阅 造林 单位 的 造林 记录 确定 林 分 
的 年 龄 (4)。 从 调查 数据 中 选择 胸径 在 5.0 cm 以 
上 生长 正常 且 未 受到 明显 损害 的 树木 ,用 于 模型 构 
建 。 本 次 研究 共 采 用 了 702 棵 樟 子 松 调查 数据 ,将 
该 数据 随机 分 为 2 部 分 ,一 部 分 是 由 17 块 标准 地 的 
501 棵 樟 子 松 组 成 的 建 模 数据 , 另 一 部 分 是 由 5 块 
标准 地 的 201 棵 樟 子 松 组 成 的 检验 数据 。 

1.3 研究 方法 

1.3.1 基础 模型 及 广义 模型 的 确定 确定 基础 模 
型 时 ,选择 12 个 常用 的 冠 幅 一 胸径 模型 ( 表 1 ) 分 别 
进行 拟 合 (最 小 二 乘法 ,0OLS ) 与 检验 ,选择 评价 指标 
决定 系数 (R ) 较 大 ,平均 绝对 误差 (MAE ) 和 均 方 根 
误差 (RMSE ) 较 小 的 模型 作为 基础 模型 。3 个 指标 
的 计算 公式 分 别 为 : 


 =1-$3(ow,-cW,)?/¥ Sow, -CW,)" 


(1) 
MAE =+3  |cw,-CW,) | (2) 
n izil 4 4 
ruse = +$ 5 (cw, -c W) (3) 
n i=lj= i 


式 中 :n 为 样本 树木 总 数量 ;m 为 标准 地 个 数 ;m 为 
第 ; 块 标准 地 内 樟 子 松 样本 株数 ,i =1,2,…,m; 
CW, All C 配 分 别 为 第 ;个 标准 地 第 7 棵 树木 冠 帆 测 
量 值 和 预测 值 (m) ;CW 为 冠 幅 测量 值 的 平均 值 (m) 。 

在 确定 基础 模型 后 ,为 使 模型 能 够 反映 冠 幅 一 
胸径 关系 随 林木 和 林 分 的 变化 ,采用 再 参数 化 的 方 
法 向 模型 中 加 入 树木 因子 和 林 分 变量 以 构建 广义 模 
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R1 冠 幅 一 胸径 关系 候选 模型 
Tab.1 Candidate models for canopy-DBH relationship 


模型 表达 式 形式 模型 表达 式 形式 
M1 CW =B +8, DBH +€ 线性 M7 CW = exp(By +8, DBH) +e Growth 
M2 CW=Bo +8, DBH +B, DBH? +E 二 次 项 M8 CW =1/( By +6, DBH) +e Hossfeld 
M3 CW =B DBH®! +€ FE PRL M9 CW = [DBH/(Bo +6, DBH) J +e Hossfeld 1 
M4 CW =B, BPH +e 复合 M10 CW =B,) (1 -exp( -BIDBH)]+e Mitscherlich 
MS CW = Byexp(B, DBH) +e 指数 Mil CW = Byexpl -Biexp( -8,DBH) J Gompertz 
M6 CW =B, +8, InDBH +€ 对 数 M12 CW =B)/(1 +8, exp( -B,DBH)] +€ Logistic 
WÈ bo B, FB, 为 模型 参数 ,s 为 误差 项 。 
表 2， 樟 子 松 调查 林 分 基本 情况 统计 
Tab.2 Basic information of the investigated forest stand of Pinus sylvestris var. mongolica 
oe 建 模 数据 由 =501 检验 数据 n, =201 
均值 标准 差 极 小 值 极 大 值 均值 标准 差 极 小 值 极 大 值 
冠 幅 /m 3.6 0.9 1.2 7.9 3.2 1.1 0.6 5.9 
胸径 /cm 16 4.4 6.4 29.1 14.5 4 6.7 27.1 
树 高 /m 9.4 2.3 3.9 16.1 9.9 1.8 6 15.5 
枝 下 高 /m 4.0 1.6 1.1 9.2 4.7 1.8 1.3 9.7 
高 径 比 0.59 0.12 0.24 1.07 0.7 0.2 0.41 1.22 
冠 长 /m 5.3 1.5 0.8 9.6 $32 1.5 2 9.6 
冠 长 率 0.58 0.1 0.25 0. 82 0.5 0.1 0.23 0. 82 
优势 木 平均 高 /m 10.7 2.1 7.33 13.61 11 231 8.2 13.2 
株数 密度 /( 株 . hm~?) 740 387 300 1 675 1 180 1 042 400 2 975 
每 公顷 断面 积 /m? 16.09 3.91 10. 04 22.29 19 8 14. 28 33.34 
相对 植 距 0.38 0.09 0.25 0.56 0.3 0.1 0.16 0.4 
平均 胸径 /cm 17.8 37 11.5 23.9 16.9 3.7 13 22.4 
算术 平均 胸径 /cm 17.6 3-7 11.2 23.7 16.6 3.9 12.7 22.2 
林 龄 /a 36.5 12.2 25 56 38.4 12.7 25 56 


型 。 在 参考 文献 [1,8 -9,21] 基 础 上 ,选择 描 
述 林木 因子 差异 和 林 分 特征 的 变量 进行 测试 。 林 木 
因子 :五 .HCB 、 树 高 和 胸径 比值 (HDR ) 、 冠 长 (CL， 
树 高 与 第 一 活 校 下 高 之 差 ) , 冠 长 率 ( CR, 冠 长 与 树 


应 模型 ,其 一 般 形式 为 : 
yi =f (B, Ui, %i) +4; (5) 
式 中 :y 为 (n; x1) 维 冠 幅 测量 值 组 成 的 因 变 量 丘 


高 比值 ) ; 林 分 变量 :优势 木 平均 高 (7H) 、 株 数 密度 
(N) 每 公 需 断面 积 (B4)、 相 对 植 中 (RS, 用 公式 
(4) 计 算 ] .平方 平均 胸径 (CQMD ) 算术 平均 胸径 
(AMD) , 林 龄 (4)。 其 中 TH 用 每 块 标准 地 内 选择 
4 ~5 棵 胸径 最 大 树木 的 树 高 平均 值 表示 ,用 以 代替 
YH SBE 。 树 木 因 子 及 林 分 变量 的 基本 情况 
见 表 2。 


V10 000/N 
R= y (4) 
式 中 :RS 为 相对 植 距 ; N 为 林 分 密度 ( 株 : hm ~*) 
TH 为 优势 木 平均 高 (m)。 


1.3.2 非 线 性 混合 效应 模型 ”本 研究 建立 的 冠 幅 
混合 模型 是 基于 标准 地 效应 的 单 水 平 非 线 性 混合 效 


量 ;%; 为 自 变量 (广义 模型 中 含有 林木 因子 和 林 分 
变量 ) ;f 为 包括 自 变 量 与 参数 的 函数 形式 ;B (p x 
1 ) 维 固定 效应 参数 向 量 ;u; 为 (q x1) 维 随机 效应 参 
数 向 量 ,wu ~N (0, D) sp 和 4 分别 为 固定 效应 参数 
和 随机 效应 参数 的 个 数 ;e; 为 误差 向 量 ,e ~N (0, 
R,) sD 为 (q xgq) 维 随机 效应 方差 - 协 方 差 矩阵 ;R， 
为 (n; xn) 维 标准 地 内 方差 - 协 方差 结构 。 

在 基础 模型 上 建立 混合 效应 模型 需要 确定 模型 
中 的 参数 效应 ,标准 地 内 方差 - 协 方差 结构 (R;) 和 
随机 效应 参数 矩阵 (DD) 。 确 定 参 数 效 应 时 ,对 所 
有 的 混合 效应 参数 组 合 进行 拟 合 ,选择 所 有 拟 合 结 
果 中 赤 池 信息 量 准 则 (AIC)、 贝 叶 斯 信息 准则 
(BIC) 和 人 负 2 倍 对 数 似 然 函数 值 ( -2LL) 最 小 的 组 


FOR 区 


合 形式 ,同时 使 用 似 然 比 检验 (ZR7) 来 避免 过 度 参 
数 化 。 在 确定 R, 时 需要 同时 考虑 观测 数据 的 异 方 
差 性 和 自 相 关 性 ,本 文 数据 未 涉及 到 自 相关 性 , 且 在 
分 析 过 程 中 未 发 现 明 显 异 方差 性 (图 1) ,因此 用 
R; =o In; 描述 标准 地 内 方差 - 协 方差 结构 。 且 最 
终 模型 中 只 有 1 个 随机 效应 参数 ,因此 ,其 随机 效应 
方差 - 协 方差 结构 为 D =o 
1.3.3 ”模型 拟 合 与 检验 ”所 有 的 数据 处 理 均 在 R 
软件 中 完成 。 混 合 效应 模型 参数 估计 使 用 R 软件 
nlme 程序 包 中 极 大 似 然 估 计 实 现 。 此 外 ,使 用 最 
小 二 乘法 拟 合 所 得 模型 简称 为 OLS 模型 。 

混合 模型 可 进行 总 体 平均 水 平 (MPA, 即 随机 
参数 设 为 0) 预测 和 具体 标准 地 水 平 (MPS ) 预测 ,其 
进行 MPA 预测 时 与 OLS 模型 相同 。 进 行 MPS 时 ， 
需 使 用 最 优 线 性 无 偏 估 计 (EBLUP) 计 算 该 标准 地 
随机 参数 值 , 其 表达 式 为 : 

ù, =DZ'(Z,DZ! +R.) ~'e, 


(6) 
RP a, 为 随机 效应 参数 的 最 优 无 偏 估计 值 ;D 为 随 
机 效应 参数 的 方差 - 协 方差 矩阵 ;R, 为 标准 地 内 方 
HE - 协 方差 结构 矩阵 ;Z 为 随机 效应 参数 设计 ( 偏 导 
数 ) 和 矩阵 ;27 为 Z, 的 转 置 矩阵 ;2 为 观测 值 减 去 用 
固定 效应 参数 计算 的 预测 值 。 

由 于 混合 模型 的 预测 精度 受 计算 随机 参数 所 用 
的 抽样 方案 影响 ,笔者 设计 了 随机 选择 1 ~ 8 棵 样本 
木 . 选 1 ~8 棵 平均 胸径 木 (12.7 ~24.0 em) . 选 1~ 
8 棵 较 小 胸径 木 (6.7 ~21.3 cm) 和 1 ~8 棵 较 大 胸 
径 木 (15.3 ~27.1 cm)4 种 抽样 方案 。 

为 了 表示 模型 的 拟 合 精度 ,在 模型 拟 合 时 评价 


指标 除 AIC、BIC 和 - 2LL 外 ,还 包括 R MAE 和 
RMSE ,模型 的 预测 检验 指标 采用 MAE 和 RMSE, 


2 结果 与 分 析 


2.1 基础 模型 

选择 的 12 个 基础 模型 在 描述 冠 幅 一 胸径 关系 
时 ,模型 的 评价 指标 相差 较 小 ( 表 3)。 其 中 ,二 次 项 
函数 (M2) 、 窜 函数 (M3) 和 单 分 子 式 ( M10) 的 拟 合 
效果 与 预测 效果 都 较 其 他 模型 略 好 ,3 个 模型 的 拟 
合 与 检验 指标 十 分 相近 ,其 中 M3 常 作 为 冠 幅 预测 
模型 ,日 可 以 描述 冠 幅 与 胸径 之 间 的 异 速生 长 关 
AE ,因此 选择 M3 为 基础 模型 ,其 表达 式 为 : 

CW, =BuDBHA + e, 


(7) 
式 中 :CW, 和 DBH, 分 别 为 第 ; 块 标准 地 内 第 7 棵 树 
的 冠 幅 (m) 和 胸径 值 (em) ;6。 和 B, 为 模型 参数 ;si 
为 模型 误差。 

为 使 模型 能 够 反映 冠 幅 和 胸径 关系 在 不 同 林木 
和 林 分 中 的 差异 ,同时 避免 因 加 入 变量 过 多 而 产生 
多 重 共 线性 和 增加 调查 成 本 ,通过 计算 比较 ( 因 情 
况 较 多 而 未 列 出 ) ,最 终 选择 了 与 模子 松 冠 幅 相关 
程度 高 且 对 模型 影响 大 的 枝 下 高 (HCB) .相对 植 距 
CRS) 和 林 龄 (4) 加 入 模型 中 ,其 表达 式 为 : 

CW, = (Bo +B, HCB, + BRS, + B;A,) DHBie, (8) 


式 中 :Bo -Bs 为 模型 参数 ;HCB; 为 第 i 块 标准 地 内 

第 7 棵 树 的 枝 下 高 (mm) RS, 为 第 i 块 标准 地 的 相对 

植 距 ;4, 为 第 i 块 标准 地 的 林 龄 (a)。 
模型 的 拟 合 结果 ( 表 4) 表明, 广义 OLS 模 型 


RI 候选 模型 拟 合 及 检验 结果 
Tab.3 Results of the fitting and test of the candidate models 


we 参数 建 模 数 据 检验 数据 

Bo Pı Bo MAE RMSE R? MAE RMSE 
M1 1.024 3 0. 163 0 0.471 4 0.607 7 0.587 2 0.620 3 0.743 5 
M2 0.824 7 0. 188 6 -0.000 8 0.470 4 0. 607 3 0.587 6 0.6170 0.738 5 
M3 0. 504 3 0.715 0 0. 469 6 0.607 5 0.587 5 0.6149 0.736 3 
M4 1.8279 1.042 8 0. 482 3 0.616 9 0.574 6 0.643 1 0.7749 
M5 1.827 9 0.041 9 0. 4823 0.616 9 0.574 6 0.643 0 0.774 8 
M6 -3.152 8 2.483 4 0.4710 0.6177 0.573 5 0.613 9 0.733 9 
M7 0.603 2 0.041 9 0. 482 3 0.616 9 0.574 6 0.643 0 0.7749 
M8 0. 442 3 -0.010 0 0. 500 4 0.636 2 0.547 5 0.676 0 0.818 7 
M9 2.957 2 0.333 5 0.469 9 0.6149 0.577 3 0.604 2 0.723 7 
M10 8.2375 0. 037 1 0.469 2 0.610 2 0.583 8 0.604 7 0.725 2 
M11 8.664 7 2.020 6 0.053 0 0.4709 0.607 7 0.587 2 0.617 2 0.739 7 
M12 7.102 7 4.383 2 0.095 5 0.471 5 0. 608 2 0.586 5 0.6175 0.740 4 
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表 4 不 同 冠 幅 一 胸径 模型 拟 合 结果 及 统计 指标 
Tab.4 Fitting results and statistic indexes of different 


canopy-DBH models 


变量 模型 (7) 模型 (8 ) 模型 (9) ”模型 (10) 
Bo 0. 504 3 0. 342 7 0. 309 5 0. 323 4 
Bi 0.7150 -0.0077 0.8882 -0. 006 9 
Bo 0.1310 0.1219 
Bs -0.001 2 -0.001 3 
Bs 0. 882 3 0.910 2 
oo" 0.001 0 0. 000 2 
o 0.3705 0.305 6 0.274 4 0.2709 
AIC 928.31 834.78 828.55 806.97 
BIC 940.96 860.07 845.42 836.49 
-2 LL 922.31 822.78 820.55 792.97 
R? 0.587 5 0.661 8 0.703 0 0.705 1 
MAE 0.469 6 0.415 5 0.386 6 0. 382 2 
RMSE 0.607 5 0.550 0 0.515 4 0.513 6 


(8) fy AIC BIC , -2 LL MAE Ail RMSE 较 基 础 OLS 
模型 (7 ) 分 别 降 低 了 10.08% .8.60% 10.79% 、 
11.51% 和 9.46% ,R° 升 高 了 12.65% 。 对 2 个 模型 
进行 方差 分 析 后 发 现 ,2 个 模型 之 间 存 在 极 显 著 差 
FF =36. 34,P <0. 001) ,可见 模型 (8) 较 模型 (7) 
更 适合 用 于 冠 幅 预测 。 
2.2 混合 效应 模型 

以 不 同 的 组 合 形 式 向 模型 (7) 和 模型 (8 ) 中 加 
入 随机 效应 参数 ,发现 2 个 模型 均 为 Bu 中 含有 标准 
地 效应 的 随机 参数 时 ,模型 的 拟 合 效果 最 好 ,2 个 混 
合 模型 表达 式 分 别 为 : 

CW, = (Bo + u;) DBH + si (9) 


CW, = (By +u; +8, HCB, + B,RS; + 
BsA;) DBH® +e, (10) 
式 中 :uw 为 随机 参数 ,其 他 变量 含义 同上 文 。 


4 (a) 模型 (9) 


冠 幅 预测 值 /m 
图 1 冠 幅 混 


沙 林 冠 幅 一 胸径 模型 


混合 模型 在 拟 合 过 程 中 都 未 出 现 明 显 的 异 方差 
性 (图 1), 因 此 不 考虑 进行 异 方差 处 理 。 模 型 
(7) ~ (10) 的 拟 合 结 果 表 明 ( 表 4) ,含有 混合 模型 
的 拟 合 精度 明显 高 于 OLS 模型 ,广义 模型 的 拟 合 效 
果 要 好 于 采用 相同 拟 合 方法 的 基础 模型 。 对 模型 
(9) 和 模型 (10 ) 进行 似 然 比 检验 ,发 现 模型 (10 ) 的 
拟 合 优 度 显 著 优 于 模型 (9) (LRT = 27. 58,P <0. 
000 1) 。 相 较 于 模型 (9) ,模型 (10) 的 尽 升 高 了 0. 
30% ,MAE 和 RMSE 分 别 降 低 了 1.14% 和 0.35% 。 
基于 模型 (9) 和 模型 (10) 的 混合 效应 , 画 出 2 个 模 
型 在 不 同 林 分 中 的 冠 幅 一 胸径 曲线 与 观测 值 散 点 分 
布 组 合 图 (图 2) ,其 中 模型 (10) 中 的 校 下 高 变量 以 
每 块 标准 地 中 枝 下 高 平均 值 代替 ” 。 如 图 2 所 示 ,2 
个 模型 拟 合 的 冠 幅 一 胸径 关系 曲线 均 穿 过 冠 幅 一 胸 
径 的 散 点 分 布 区 域 ,多 数 集中 于 中 心 区 域 ,少量 靠近 
较 边 缘 区 域 。 说 明 混 合 模 型 具有 较 高 的 精度 ,可 以 
很 好 地 模拟 不 同 樟 子 松林 分 中 冠 幅 一 胸径 的 关系 。 
2.3 ”模型 预测 与 抽样 对 比 

由 模型 预测 精度 和 抽样 方案 关系 可 知 (图 3):4 
种 抽样 方案 中 ,选择 平均 木 时 模型 预测 精度 最 高 , 且 
MAE 和 RMSE 都 随 样本 木 数 量 的 增加 而 出 现 减 小 
趋势 。 对 于 广义 混合 模型 ,其 预测 精度 整体 变化 较 
小 , 且 MAE 和 RMSE 略 低 于 基础 混合 模型 。 对 于 基 
础 混合 模型 :选择 2 棵 平均 木 时 ,其 MAE 和 RMSE 
分 别 较 其 MPA 预测 降低 了 28. 1% 和 20.9% ,之 后 
随 样 本 木 数 量 增加 降低 速度 较 小 ;用 同样 的 方法 发 
现 其 余 3 种 方案 中 样本 木 数量 大 于 4 棵 时 ,MAE 和 
RMSE 的 降低 速度 变 绥 。 此 外 ,本 研究 中 模型 (9 ) 采 
用 随机 选择 、 选 小 树 和 选 大 树 3 种 方案 时 ,其 MPS 
预测 效果 不 如 模型 (8 ) 和 模型 (10) 的 MPA 预测 。 


4 (b) 模型 (10) 
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Fig.1 Residual distribution of canopy mixed effect models 
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(a) 模型 (9) (b) 模型 (10) 


胸径 /cm 胸径 /cm 
图 2 和 冠 幅 一 胸径 观测 值 与 不 同 标准 地 拟 合 曲线 对 照 


Fig.2 Comparison between the canopy-DBH observed values and the sample-plot fitted curves 


一 模型 (9) -一 模型 (7) ---M9 一 模型 10) 一 模型 (8) ---M10 
0.7 | (a) 随机 选择 r (b) 选 平 均 木 - (© 选 小 树 (d) 选 大 树 
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注 :M9 和 M10 分 别 表示 模型 (9) 和 (10) 的 MPA 预测 。 
图 3 模型 预测 精度 与 样本 木 数 量 关 系 


Fig.3 Relationship between the prediction accuracy and the number of sampling tree 
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图 4 不 同 条 件 下 冠 幅 一 胸径 曲线 模拟 


Fig.4 Canopy-DBH simulation curves under different conditions 
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韩 艳 刚 等 : 樟 子 松 人 工 固 沙 林 冠 幅 一 胸径 模型 


2.4 ， 樟 子 松 人 工 林 冠 幅 模拟 

使 用 模型 (10) 的 MPA 预测 功能 对 不 同 林 分 和 
个 体 的 樟 子 松 冠 幅 进行 模拟 ,模拟 分 析 DBH 之 外 的 
某 个 增加 变量 (HCB、RS 和 4) 对 冠 幅 预测 的 影响 
时 , 除 目 标 变量 外 ,其 余 增 加 变量 均 采 用 调查 数据 的 
平均 值 代 替 。 由 图 4 可知, 在 使 用 模型 (10 ) 进行 冠 
幅 预 测 时 ,加 入 模型 中 的 林 分 变量 和 林木 因子 ,对 冠 
幅 预测 的 影响 与 其 参数 估计 值 的 符号 存在 着 明显 的 
对 应 关系 , 冠 幅 随 着 HCB 的 增 大 而 明显 减 小 ; 冠 幅 
与 RS 的 关系 与 HCB 相反 ,表现 为 冠 幅 随 RS 的 增 
大 而 增 大 ;从 樟 子 松 由 中 龄 林 到 过 熟 林 的 生长 变化 
过 程 来 看 , 冠 幅 随 4 的 增加 而 呈 减 小 趋势 。 


3 讨论 


3.1 冠 幅 模型 比较 

本 研究 以 章 古 台地 区 中 龄 林 到 过 熟 林 中 不 同 立 
地 和 不 同 密度 的 樟 子 松 人 工 纯 林 为 例 , 从 12 个 常用 
冠 幅 一 胸径 中 选 出 了 拟 合 与 检验 效果 较 好 的 震 函 数 
形式 作为 基础 模型 ,并 从 多 个 树木 因子 和 林 分 变量 
中 筛选 出 了 枝 下 高 (HCB ) 、 相 对 植 距 (RS) 和 林 龄 
(4) 作 为 模型 的 增加 变量 。 在 寡 函 数 模型 基础 上 ， 
检验 了 增加 预测 变量 和 建立 混合 效应 模型 对 模型 精 
度 的 影响 ,结果 表明 ,建立 混合 效应 模型 和 广义 模型 
都 可 以 提高 模型 的 精度 ,2 种 方法 结合 后 的 广义 混 
合 模型 (10 ) 的 精度 最 高 ,但 与 含有 胸径 单个 自 变 量 
的 简单 混合 模型 (9) 之 间 差 异 不 大 。 

通过 比较 不 同 模型 在 预测 检验 时 的 MAE 和 
RMSE, 发现 混合 模型 MPS 预测 精度 高 于 相应 OLS 模 
型 预测 。 模 型 (10) 的 预测 精度 略 高 于 模型 (9 ) ,但 两 
者 差异 不 大 ,类似 的 情况 也 出 现在 了 一 些 树 高 一 胸径 
模型 的 研究 中 3 ,这 可 能 是 由 于 模型 中 的 随机 参 
数 已 经 能 够 反映 不 同 林 分 和 树木 间 的 差异 ,此 时 
加 入 树木 因子 和 林 分 变量 虽 可 提高 模型 的 精度 *， 
但 可 能 作用 不 大 。 混 合 模型 的 MPA 预测 与 OLS 
模型 预测 相似 ,其 精度 主要 受 自 变 量 影响 。 由 于 广义 
模型 中 的 林木 因子 和 林 分 变量 能 够 反映 林木 和 林 分 
的 差异 ,因而 其 精度 较 简单 ,模型 精度 高 。 

研究 表明 ,不 同 树种 适用 于 不 同 的 抽样 设计 方 
法 。 本 研究 从 设计 的 4 种 抽样 方案 中 发 现 ,广义 
混合 模型 和 基础 模型 的 精度 都 随 样本 木 数 量 的 增加 
而 升 高 , 且 使 用 平均 木 时 模型 的 预测 精度 最 高 。2 
个 混合 模型 在 使 用 2 棵 平均 木 计算 随机 参数 时 , 预 


测 精度 MAE 和 RMSE 的 降低 速率 较 快 , 增 加 平均 木 
的 数量 可 以 提高 模型 精度 ,但 是 变化 速率 较 慢 , 且 需 
增加 调查 成 本 。 因 此 ,推荐 使 用 2 棵 平均 木 计算 混 
合 模型 随机 参数 值 。 

林 分 结构 与 立地 差异 较 小 时 ,基础 OLS 模型 可 
以 实现 相似 林 分 的 准确 预测 。 林 分 间 差 异 较 大 
时 ,混合 模型 和 广义 OLS 模型 都 能 够 反映 林 分 和 树 
木 差 异 对 冠 幅 的 影响 。 若 无 平均 木 冠 幅 测量 值 作为 
先 验 信 息 ,利用 其 他 方案 选择 较 少 样本 木 计算 基础 
混合 模型 的 随机 参数 值 ,会 导致 其 预测 精度 低 于 广 
义 模 型 , 旦 广义 混合 模型 MPA 、MPS 及 广义 OLS 模 
型 间 差 异 较 小 ,此 时 使 用 广义 OLS 模型 (8 ) 进行 冠 
幅 预 测 更 为 简便 且 精 度 较 高 ,类似 结果 也 出 现在 了 
Gomez-Garcia 等 1 研究 中 
3.2 ” 冠 幅 与 树木 因子 及 林 分 变量 关系 

本 研究 所 加 变量 (HCB RS 和 4) 均 与 模子 松 冠 
幅 存在 极 显 著 相关 性 。 其 中 ,HCB 是 反映 树冠 衰退 
和 树 高 生长 动态 的 重要 因子 之 一 ,对 于 喜光 树种 樟 
子 松 来 说 , 树 高 的 生长 会 导致 下 层 枝叶 接受 光照 较 
少 而 死亡 ,这 种 现象 在 密度 大 的 林 分 中 更 明显 ; 当 树 
冠 发 生 衰退 时 ,树木 活力 下 降 ,一 般 下 层 枝 条 先 死 
亡 , 此 时 HCB 升 高 , 且 常 出 现 冠 幅 减 小 。 因 此 ， 
常 使 用 HCB 作为 冠 幅 的 预测 变量 , AL oe BY BE 
HCB 升 高 而 变 小 。 对 于 林 分 来 说 ,RS 是 一 个 描述 
林 分 密度 的 综合 指标 , 可 能 是 该 指标 在 计算 时 包含 
了 描述 林 分 竞争 的 N 和 代表 立地 质量 的 7 五 ,本 文 
在 将 THN 和 BA 加 入 模型 中 时 ,其 效果 均 没 有 加 入 
RS 时 好 。 在 林 分 中 ,RS 的 增 大 意味 着 树木 可 占有 
的 相对 生长 空间 增 大 ,有 利于 树冠 横向 生长 , 导 
致 冠 幅 随 RS 的 增 大 而 增 大 。 除 ACB 和 RS 之 外 ， 
还 加 入 了 4 来 表示 林 分 特征 。 在 造林 资料 保存 完 
整 的 单位 中 ,人 工 纯 林 的 林 龄 是 最 易 获取 的 林 分 指 
ER ,根据 林 龄 不 仅 可 以 对 树 高 ,胸径 和 和 冠 幅 进 行 预 
测 , 还 可 以 整体 把 握 林 分 动态 。 由 于 林 龄 与 其 
他 多 个 因子 (特别 是 TH) 之 间 存 在 较 强 的 相关 性 ， 
因而 常 在 建立 广义 模型 时 被 舍 去 2 ,但 本 文 结 果 表 
明 ,在 含有 树木 因子 HCB 和 林 分 变量 RS 的 模型 中 ， 
加 入 林 龄 ,对 模型 精度 的 提高 要 高 于 其 他 变量 , 故 将 
林 龄 作为 冠 幅 预测 变量 之 一 。 对 于 树木 来 说 , 冠 幅 
的 生长 通常 是 一 个 随 着 年 龄 增长 而 增 大 的 过 程 。 本 
文 在 检验 林 龄 的 影响 时 ,结果 却 是 冠 幅 在 中 龄 林 到 
过 熟 林 的 过 程 中 会 逐渐 减 小 ,这 可 能 是 由 于 樟 子 松 
的 冠 幅 生 长 主要 集中 于 幼 龄 林 和 中 龄 林 阶 段 ,35 a 
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左右 会 出 现 提早 衰弱 ”' ,使 樟 子 松 人 工 林 在 中 龄 林 
之 后 整体 表现 出 冠 幅 随 林 龄 增 大 而 减 小 的 现象 。 但 
是 关于 和 冠 幅 随 着 林 龄 增长 而 减 小 的 内 在 原因 ,还 需 
进一步 的 研究 验证 。 

4 结 
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本 文通 过 选择 合理 的 沙 地 樟 子 松 人 工 林 标 准 地 
进行 调查 ,利用 调查 数据 构建 了 章 古 台地 区 的 樟 子 
松 人 工 林 冠 幅 预 测 模型 ,得 出 了 如 下 结论 : 

(1) 老 函 数 模型 可 以 较 好 地 描述 章 古 台 樟 子 松 
人 工 林 的 冠 幅 一 胸径 关系 ,适合 作为 预测 冠 幅 的 基 
础 模型 。 

(2) 混合 效应 模型 和 广义 模型 可 提高 冠 幅 模 型 
精度 ,但 是 广义 混合 模型 .基础 混合 模型 及 广义 OLS 
模型 三 者 在 用 于 冠 幅 预测 时 差异 不 大 , 故 推荐 使 用 
所 需 变量 少 的 基础 混合 模型 或 方法 简便 的 广义 OLS 
模型 进行 冠 幅 预测 。 对 于 基础 混合 模型 ,建议 每 块 
标准 地 中 选择 2 棵 平均 木 计算 随机 参数 值 。 

(3) 所 考虑 的 树木 因子 和 林 分 变量 中 , 校 下 高 、 
相对 植 距 和 林 龄 对 冠 幅 预测 模型 中 的 影响 最 大 ,可 
以 显著 提高 模型 精度 ,而 且 冠 幅 一 胸径 曲线 与 这 些 
变量 之 间 的 关系 可 以 通过 广义 模型 模拟 实现 。 
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Canopy-DBH Models Forsand-fixing Plantation of Pinus sylvestris var. mongolica 


HAN Yan-gang', LEI Ze-yong', ZHAO Guo-jun’, ZHOU Yan-ping', XU Chang' 
College of Environmental Science and Engineering ,Liaoning Technical University , Fuxin 123000 , Liaoning , China; 


(1. 
2. Liaoning Provincelnstitutefor Sand-fixation and Afforestation Research , Fuxin 123000 , Liaoning , China ) 


Abstract; In this study, the data consisted of 702 individuals of Pinus sylvestris var. mongolicain 22 temporary 
sample plots of pure plantation in Zhanggutai of Liaoning Province , northeast China, were used to develop the cano- 
py-DBH models, i. e. the basic model ( fitted by ordinary least squares (OLS) ) , generalized model (fitted by 
OLS) ‚nonlinear mixed-effect basic model and nonlinear effect-mixed generalized model. The goodness-of-fits and 
prediction accuracy of the four models were compared. For the effect-mixed model, four sampling strategies ,i. e. , 
random sampling , large-DBH tree sampling, small-DBH tree sampling and medium-DBH tree sampling, were de- 
signed to calculate the random parameters. The effects of individual factor and stand level variables on the canopy- 
DBH relationship were simulated. Model evaluation indices included the determination coefficient ( R°) ,mean abso- 
lute error (MAE) and root mean square error (RMSE). Results showed that the height to live crown base (HCB) , 
relative spacing index (RS) and stand age (A) were the dominant factors in CW-DBH models. The goodness-of-fits 
of effect-mixed CW-DBH models (R? ,MAE and RMSE) were 0.703 0,0. 386 6 and 0.515 4,and those of effect- 
mixed generalized model were 0.705 1 ,0.382 2 and 0.513 6,respectively , which were better than those of the OLS 
models. However, the difference in goodness-of-fit between the effect-mixed basic and generalized models was not 
significant (below 1% ). Canopy decreased with the increase of A and height of HCB, but increased with the in- 
crease of RS. 
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